Gundlach oscillations and Coulomb blockade of Co nano-islands on MgO/Mo(100) 
investigated by scanning tunneling spectroscopy at 300 K 
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Ultrathin MgO films on Mo(100) with a thickness up to 12 ML are studied by scanning tunneling 
microscopy and spectroscopy at room temperature. The spatial variation of the work function within 
the MgO film is mapped by field emission resonance states (Gundlach oscillations) using dz/dU 
spectroscopy. We found circular spots with significantly reduced work function (A$ = 0.6 eV), which 
are assigned to charged defects within the MgO film. On top of the MgO films, small Co cluster are 
deposited with an average contact area of Ac — 4nm 2 . These islands exhibit Coulomb oscillations 
in dl /dU -spectra at room temperature. Good agreement with orthodox theory is achieved showing 
variations of the background charge Qo for islands at different positions, which are in accordance 
with the work function differences determined by the Gundlach oscillations. 

PACS numbers: 73.23.Hk, 73.40.Gk, 68.55.aj, 48.47.Gh 



I. INTRODUCTION 



Magnesium oxide (MgO) is a preferred insulator for 
magnetic tunnel junctions p]-[9| and is extensively used 
as a template for the microscopic study of catalytic re- 
actions pjj-[21[. Both fields rely on the good epitaxial 
quality of MgO films, while the former was additionally 
boosted by the theoretical insight of a symmetry induced 
strong spin selectivity of MgO(OOl) [22]. MgO might also 
be an excellent template for nanoelectronic studies by 
scanning probe microscopy similar to the ones that have 
been performed recently on NaCl with resp ect to charge 
manipulation and bond formation |23l-[25| or on AI2O3 
and CuN with respect to the determination of magnetic 
properties of individual atoms on a substrate [26|-[28j . 
Thin films of MgO, exhibiting a wide band gap and a sim- 
ple rock salt structure, grow epitaxially on different metal 
substrates as Ag(001) QBfifil, Fe(001) 0-0 and Mo(001) 
[l7l-[2lj. Here, we choose Mo, since it allows high anneal- 
ing temperatures exceeding 1000 K, which might foster 
an improved MgO film quality. The remaining major de- 
fects are color centers exhibiting defect states within the 
band gap with energies depending on the defect position 
jl3-KL6j . Work function differences between different layer 
thicknesses of MgO have been determined by field emis- 
sion resonance (FER) spectroscopy and photon mapping 
to be up to 650 meV |2fl|. 

In this work, we employ scanning tunneling microscopy 
(STM) and spectroscopy (STS) in order to investigate 
pure MgO films grown epitaxially on Mo(100) as well as 
Co nano-islands deposited on top of the MgO. We observe 
an increase of the MgO band gap Eq&p with increasing 
film thickness which reached EQ ap — 7.1 eV at a thick- 
ness of 11 monolayers (ML). Mapping the work function 
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y) by the energetic shift of the first Gundlach os- 
cillation reveals a strong spatial variation up to about 
A<£ = 0.6 eV. This shift could be correlated with a peak 
within the gap known to be caused by charged color cen- 
ters [14{ . This explains the work function shifts straight- 
forwardly as caused by color centers located within the 
topmost two MgO layers. 

Probing Co islands on top of the MgO film, we observe 
Coulomb oscillations by dl/dU spectroscopy with ener- 
gies in excellent agreement with orthodox theory (29j . 
A varying Coulomb gap around U = 0V indicates a 
fluctuating background charge, which could also be at- 
tributed to the charged color centers close to the surface. 
The width of the Coulomb peaks of 0.3 — 0.4 eV is much 
larger than the energy resolution of the experiment. This 
is probably due to a large coupling of the cluster charge 
to the phonons of the ionic insulator [30]. Notice that 
Coulomb oscillations of metal clusters on insulating thin 
films have rarely been observed at room temperature so 
far @,[3i|, 11. 



II. EXPERIMENT 

The experiments were carried out using an ultrahigh 
vacuum system (base pressure p = 8-10 -11 mbar) consist- 
ing of two separately pumped chambers. The main cham- 
ber is equipped with a 4-grid low-energy electron diffrac- 
tion (LEED) optics, an electron beam heating stage and 
a modified Omicron STM for measurements at room 
temperature. A side chamber is equipped with the fa- 
cilities for MgO growth, namely an electron beam in- 
duced Mg evaporator, a leak valve for oxygen and a 
heating stage positioned directly in front of the evapora- 
tor. Prior to MgO preparation, the Mo(100) crystal was 
cleaned by cycles of annealing in an oxygen atmosphere 
(Po 2 — 5 • 10 _T mbar) at 1400 K followed by subsequent 
flashing to 2300 K. MgO films were prepared by molec- 
ular beam epitaxy of magnesium in an O2 environment 
at a partial pressure of po 2 = 1.5 • 10 _7 mbar keeping 
the substrate at room temperature. The deposition rate 



of 0.5ML/min has been controlled by means of a quartz 
microbalance operating at 10 MHz. After MgO deposi- 
tion, the samples were annealed at HOOK for lOmin. We 
checked the purity of the substrate and the MgO films 
by means of Auger electron spectroscopy and LEED. All 
STM and STS measurements were carried out with an 
electrochemically etched tungsten tip with the bias ap- 
plied to the sample. STS data (dl/dU, dz/dU) were ac- 
quired by lock-in-technique applying a modulation volt- 
age with amplitude V^'^mod- dl/dU curves are recorded 
with feedback off after stabilizing the tip at voltage U s t&b 
and current 7 s tab- 



III. RESULTS AND DISCUSSION 

Figure Q] a) shows an STM image of an MgO film 
on Mo(100) with a thickness of 7-8 ML. The nearly 
completely covered surface exhibits 1.3 nm deep grooves 
(black lines) which mainly separate differently tilted 
MgO facets. A Moire pattern caused by the 5.4% 
lattice-misfit of Mo and MgO is clearly visible all 
over the film. Previous LEED studies indicated this 
superstructure to exist up to a MgO thickness of 12 ML 
[21] , but it has never been observed for such thick 
films directly by STM. In addition, we observed large 
variations in the periodicity of the Moire pattern ranging 
from 3nm to 8.3 nm compared to the calculated value 
of 5.35 nm. Further different MgO facets show tilting 
angles with a deviation of up to 20° from the Mo [001] 
direction as indicated by the lines in Figure [1] a). This 
may be caused by relaxation effects during cooling down 
the sample from HOOK. Figure [1] b) displays a zoom 
into the film highlighting that the grooves sometimes 
appear at step edges originating from screw dislocations 
within the film. 

Figure [1] c) shows another MgO film, nominally pre- 
pared under the same conditions as the film shown in 
a) and b). Here, we do not find a smooth film, but 
large valley areas covered by only 1-2 ML MgO and, 
in between, higher MgO islands with thicknesses up 
to 11-12 ML. The thickness of the MgO within the 
valleys is consistently determined by observation of step 
heights within the valleys up to 2 ML and comparison 
of the total MgO film topography with the nominally 
deposited amount of Mg. The origin of this different 
growth type is not known, but may be caused by a 
small amount of molybdenum oxide remaining on the 
substrate prior to MgO preparation, which might alter 
the MgO nucleation. However, we could not detect 
the corresponding oxygen on the cleaned substrate by 
Auger electron spectroscopy. The Moire superstructure 
induced by misfit can still be observed on MgO islands 
with a thickness of 11 ML (see inset of figure (He). 

The band gap of the MgO film is measured by STS 
as shown in figure [1] d) . The edge of the band gap is 
determined by the point where positive dl/dU intensity 
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FIG. 1: (Color online) a) (200 x 200) nm 2 STM image of 
7-8 ML MgO film on Mo(100) (U = 3 V, I = 0.5 nA); the 
typical Moire pattern caused by the lattice mismatch of Mo 
and MgO is visible; the marked crystallographic directions 
are determined by LEED. The lines emphasize the deviation 
of the Moire pattern from the Mo [001] direction within dif- 
ferent MgO facets, b) Zoom into the MgO film as marked 
by a dashed line in image a); the Moire pattern and a screw 
dislocation indicated by an arrow is visible, c) MgO sam- 
ple consisting of 11-12 ML thick islands and valleys in be- 
tween consisting of 1-2 ML MgO. The inset shows the Moire 
pattern still appearing on MgO islands of 11 ML. d) Spec- 
tra of dI/dU(U) measured on MgO areas of 2 ML, 7 ML and 
11 ML thickness as indicated (U s tab = 3 V, 7 s tab = 0.45 nA, 
Umod = 40 mV); the deduced band gap is marked. 
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FIG. 2: (Color online) a) (40 x 40) nm 2 STM image of 7-8 ML thick MgO film (U = 3.5 V, I = 0.5 nA). b) dz/dU-map recorded 
at U = 3.5 V of the same area as shown in a) (/ = 0.5 nA, U m od = 60 mV); only the grooves in between the islands exhibit 
spectroscopic contrast, c) dz/dU map recorded at U = 4.0 V of the same sample area as a) and b) (i~ = 0.5 nA, U m od = 60 mV); 
additional bright spots appear on parts of the MgO film, d) Sketch of the potential scheme within the tunneling barrier 
explaining the field resonance states; Ef, Ey and $ marked at tip and sample are the Fermi levels, vacuum levels and work 
functions, respectively; the black line in the vacuum region is the vacuum level influenced by applied electric field and the 
image potential of tip and sample (white line without image potential); FER states are marked by n = 1, n = 2, n = 3 and 
n — 4. e) dz/dU -spectra taken at the positions marked in c) (E/gtab = 3.5 V, / sta b = 0.5 nA, U mo d = 100 mV); the first two FER 
states n = 1 and n — 2 are labeled for curve 1; the n — 1 state shifts by about AE = 620 me V between curve 1 and curve 2. 
Inset: three FER peaks measured on a dark region of c). f) dI/dU(U) -spectra recorded on two of the positions marked in c) 
(Ustah = 4.0 V, /stab = 0.5 nA, U m od = 26 mV) ; the apparent band gap exhibits a peak at U = —2.7 V only in curve 1, which 
indicates the defect state. 



is larger than the noise level. We checked that this 
method is sufficiently reliable by recording dl/dU curves 
at different stabilization parameters. An increase of the 
band gap from 5 eV on 2 ML MgO (measured within the 
valley area) towards 7.1 eV on 11 ML MgO (the high 
islands) can be deduced approaching the bulk value of 
7.8 eV [34| . Notice that the gap grows asymetrically 
indicating that electron affinity of MgO is less affected 
by film thickness than ionization energy. 
Electronic properties of the smooth film are further 
investigated by STS. Field emission resonance states 
(FER), i.e. Gundlach oscillations, are routinely observed 
as shown, e.g., in the inset of figure [2] e) exhibiting 
peaks at voltages of 4.5 V, 5.7 V and 6.4 V. Such FER 
states are known to appear if the bias voltage exceeds 
the work function of the sample. The origin of these 



states is sketched in figure [2] d). They are tunneling 
resonances caused by the constructive interference of 
electron waves, which are reflected back and forth 
between the potential step at the sample surface and 
the potential wall given by the vacuum level. Due to 
the high electric field at the required voltages to probe 
FER states by dI/dU(U), tip changes are frequent. To 
avoid that, we measured dz /(/[/-spectra with a closed 
feed-back loop keeping the tunneling current constant. 
This retracts the tip automatically at higher voltage, 
thereby limiting the electric field. We used lock-in 
technique to determine the dz/dU signal. Figure [2] 
a) shows a constant current image of a 7 ML thick 
MgO film. The corresponding dz/dU image taken at 
3.5 V (figure [2] b) shows spectroscopic contrast only at 
the deep grooves (bright lines). Figure [2] c) displays 
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the dz/dU signal of the same area taken at a bias 
voltage of 4 V, where bright spots randomly distributed 
over the MgO film appear additionally to the fading 
lines originating from the grooves. The reason for the 
spectroscopic contrast can be understood by recording 
single point dz/dU spectra on bright and dark areas. 
Figure [2] e) shows three spectra mapping the transition 
from a bright spot to the dark surrounding. The first 
FER state shifts continuously from U(n = 1) = 4 V 
(spectrum 1), measured in the center of the bright spot, 
over U{n = 1) = 4.3 V (spectrum 3) to U(n = 1) = 4.6 V 
(spectrum 2), measured on the dark surrounding 10 nm 
away from the center of the bright spot. The maximum 
shift in peak energy of AEq = 600 meV indicates a large 
change in the MgO work function. Note, that the energy 
shift of the first FER peak gives only an approximation 
of the work function difference because of the additional 
image potential, which reduces the FER energies by a 
value which is stronger at lower FER energy. Therefore, 
the peak-shift is only a lower limit of the real change of 
the work function. 

Furthermore, dl/dU spectroscopy of the band gap 
region is measured on bright spots and dark parts of the 
dz/dU image, as exemplified in figure [2] f). The dl/dU 
curves show an additional peak at about —2.7V only on 
the bright spots. Referring to theoretical calculations 
this occupied state can be attributed either to charged 
F + color centers or to MgO surface defects at 3-fold 
coordinated surface sites [33]. The defect state around 
—2.7V has been previously observed at MgO-islands on 
Ag(001) 0, [l6j, but in contrast to one of these works, 
we do not observe an empty state at «+l V. This might 
be due to a shift of the empty state towards the conduc- 
tion band of the MgO within the MgO/Mo(100) system. 
It might be surprising that the defects mapped by FER 
spectroscopy appear relatively large, i.e. about 3 — 4nm 
in diameter. This size is definitely larger than the defect 
state mapped at —2.7V in previous low temperature 
STM measurements [14|. The major difference between 
the two methods is that we detect the electric potential 
of the charge and not the state of the defect. Estimating 
the potential of a point charge about 1.5 nm away from 
the center still results in a potential shift of 200 meV. 
Together with our energy resolution of AE = 170 meV 
in Fig. [2]c), this already justifies the observed size of the 
charge within FER images. However, the field emission 
process is more subtle depending on details of the 3D 
potential between tip and substrate j35|, [36| and more 
detailed calculations are required for a more quantitative 
analysis of the apparent size of a charged defect in FER 
images. Of course charged defects of the size of a few 
atoms cannot be excluded. 

The depth of the singly positive charge with respect to 
the MgO surface (da) can be estimated by the Coulomb 
potential of a point charge (e • U c , e: elementary charge) 
leading to dd = e/(4:7rseoU c ) = 0.25 nm (so- dielectric 
constant of vacum), if the bulk value of the dielectric 
constant of MgO e = 9.8 [37| is used and dd — 0.5 nm if 



half the value of e is used as a lower estimate considering 
the presence of the surface. This indicates that the 
corresponding defect must be located very close to the 
MgO surface, probably in the subsurface layer of the 
insulating film. Notice, that previous FER spectroscopy 
and photon maps of MgO films on Mo (001) showed 
energy shifts in the same range (650 meV), but between 
thick and thin MgO-islands [20[ and not as spatial 
fluctuations within a smooth film. 
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FIG. 3: (Color online) a) Schematic drawing of the tun- 
neling junction consisting of tip/vacuum/Co- island and Co- 
island/MgO/Mo(100) together with equivalent circuit dia- 
gram, b) STM image (25 x 20) nm 2 of Co nano-islands 
on MgO (U = 3.5 V, I = 0.5 nA). c)-f) dl /dU(U) -spectra 
recorded on Co islands with different sizes Aco as indicated 
(Z/stab = 3.5 V, /stab = 0.7 nA, U mo d = 40 mV); the Coulomb 
gap AUo and the charging voltage AU are marked in c); spec- 
trum f) is measured on a Co island on 10 ML MgO, spectra 
c)-e) are recorded on an island of 11 ML MgO. 



Next, we prepared Co nano-islands on top of the MgO 
film by molecular beam epitaxy at room temperature. 
As illustrated in figure [3] b), Co forms compact islands 
exhibiting an average area of Aq — 4 nm 2 and a height 
of 2-4 ML. The Co nano-islands are well separated from 
each other. This is an ideal system for Coulomb block- 
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ade, since it consists of two tunneling junctions given by 
W-tip/ vacuum/Co-island and Co-island/MgO/Mo(100). 
The equivalent electrical circuit is drawn in figure [3] a) 
consisting of two capacitors (Ci, C2) and two resistors 
(jRi, R2) in parallel. Single electron tunneling occurs, if 
the charging energy per electron is much larger than the 
thermal energy, i.e. 



Ci+C 2 



» kT. 



(1) 



Thus, the capacity must be as small as a few atto- 
Farad in order to observe Coulomb staircases in the I(U)- 
signal, respectively, peaks in the differential conductivity 
dI/dU(U) at T = 300 K. Moreover, sequential electron 
tunneling requires resistances of the two junctions higher 
than the resistance quantum h/2e 2 = 12.9 kO, which, 
however, is easily achieved within our setup. Figure [3] 
c)-f) show dI/dU(U) -spectra for several Co islands with 
different sizes. Within the band gap of MgO (from — 5 V 
to +2V), several peaks are visible indicating Coulomb 
staircases. The central Coulomb gap AUo and the dis- 
tance between adjacent peaks at higher energy AU is 
marked in Fig. [3t). The resulting values of AUo are 
given in Fig. [3h)-f ) together with the contact area Aq 
of the corresponding island. 

To compare the measured AU with orthodox theory [29L 
we calculate the expected AU* = e/(Ci + C2) |3ll l38j. 
The capacities C\ and C2 of the two tunneling junctions 
are regarded as two plate capacitors with the area of the 
Co nano-island Aq q resulting in C = S£oAc /d. Thus, 
C2 can be determined using the measured thickness of 
the MgO-film as d and the dielectric constant of MgO 
^MgO = 9.8 [37]. For calculation of Gi, the absolute 
tip-sample distance zo must be known. This number 
is difficult to determine directly by STM. We estimate 
zo by extrapolating I(z) curves to the distance, where 
the conductance would reach the conductance quantum 
Go = 2e 2 /h. This point is regarded as the point of me- 
chanical contact corresponding to d = nm [39ll40j|. This 
results in a tip-surface distance zq at I s tab and U s t&b of: 



J- in 

2k 



1 L 



stab 



Go t^stab J 



(2) 



The deca y c onstant n for a planar tunneling junction 
is given by [41] : 



$Co - |eC4tab|)- 



(3) 



We used the work function of the densely packed 
W(110) surface <$> w = 5.25 eV for the W tip and a Co 
work function of &co — 5eV [42|. This leads to zo ~ 6.7 
A, which is taken as d for G2. Table [I] shows the cal- 
culated values Ci, C2 and AU* in comparison with the 



TABLE I: Table of calculated values for the serial capacities 
Ci, C2 and the Coulomb peak distances AU* compared with 
the measured values AU as taken from figure [3] In addition, 
the lateral size of the Co island, the determined background 
charge Qo, the corresponding local change of Co work function 
A<E>co and the average full width at half maximum of the 
Coulomb peaks AE P is given. 



size 



Co 



AU AU* \Q \/e A$ Co AE P 



3.2 nm^ 41 zF 120 zF 1.03 V 1.00 V 0.19 385 meV 395 meV 
3.7 nm 2 47 zF 139 zF 0.90 V 0.86 V 0.26 453 meV 380 meV 



4.2 nm 2 53 zF 157 zF 0.84 V 0.76 V 



meV 385 meV 



5.3 nm 2 67 zF 219 zF 0.50 V 0.56 V 0.34 408 meV 290 me V 



measured AU for different island sizes. Excellent agree- 
ment between AU and AU* is found without any fit pa- 
rameter. 

The Coulomb gap AUo depends, in addition, on the back- 
ground charge Qo and is given by [29|, [38| : 



AC/o = (| - |Qo|) 



Ci + C 2 



(4) 



Thereby, Qo can be related to the local work functions 
of the three metals according to [38] : 



Qo 



[Ci($w - $Co) - G 2 ($ Mo - $Co)] • (5) 



The determined values for Qo are also summarized 
in table HI They vary from cluster to cluster in the 
range of up to AQo = 0.34e. The corresponding spatial 
variation of the Co work function can be calculated 
as A<3>co = eQo/(C2 — Ci), if one assumes that $Mo 
is spatially constant and that the tip did not change 
between individual measurements. The former is a 
good assumption, if the background charge is located 
close to the MgO-Co interface. The resulting values of 
A<E>co are added to table H Interestingly, the maximum 
change amounts to A&co — 450 meV, which is rather 
close to the work function variation of the pure MgO 
film determined by FER spectroscopy. Thus, we assume 
that the variation of AUo also results from positively 
charged color centers located close to the MgO surface. 
The slightly lower value of 450 meV probed by the Co 
clusters with respect to 600 meV measured on the pure 
MgO might be caused by screening of the color centers 
via Co islands, by a lateral distance between Co clusters 
and color centers and/or by averaging over the lateral 
size of the cluster. 

Finally, we will discuss the full width at half maximum 
(FWHM) of the Coulomb peaks AE? determined after 
subtracting a linear background from the spectra in 
figure [3b) -f). The values of AE? averaged over the peaks 
of one curve are displayed in table HI The peak width 
of about 300 — 400 meV is much larger than the energy 
resolution of the STS experiment, which is limited 
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by the thermal energy (T = 300 K) and the applied 
modulation voltage (Z7 mo d = 40 mV) and given by 
AE T ~ ^(3.3 • kT) 2 + (2.5 • eU mod ) 2 = 130 meV 0. 
Notice that this formula is in agreement with experiment 
down to AEp = 0.1 me V accordingly T = 0.3 K [44(. In 
order to cross-check, we varied the modulation voltage 
from 30 mV to 70 mV on a single Co island observing 
changes in the FWHM of less than 20 mV. Life time 
broadening of the peaks AEl > h/2RC is orders of 
magnitudes smaller (in the range of tenth of fieV) and 
can also be neglected. Thus, the charge within the 
cluster is well equilibrated before leaving the cluster. 
Similar FWHMs (270 mV) have been observed previ- 
ously for defect states of positively charged CI vacancies 
in NaCl films even at T = 5 K [30]. They are attributed 
to an excitation of optical phonons in the NaCl lattice. 
Following this argument, we speculate that AEp is 
caused by the strong coupling of the tunneling electrons 
to optical phonons of the MgO-lattice. An estimate of 
the resulting FWHM is given by AE P = V8 In 2 • S • hu 
with S being the Huang- Rhys factor and uj being the 
relevant phonon frequency [30[. Using a typical optical 
phonon energy of MgO (Hlj ~ 50meV) [45] and the 
Huang-Rhys factor S = 39 determined for color centers 
in bulk MgO [Ifj, one gets AE P ~ 700 meV, which 
is larger by a factor of two than the measured value. 
The larger value is plausible regarding the stronger 
electron-phonon coupling within the color center than 
the distant coupling taking place between electrons 
within the Co island and phonons within the MgO. 
Alternatively, the observed width could also be due to 
charge fluctuations within the neighboring MgO, which 
are fast with respect to the measurement time, but slow 
on the time scale of a single electron tunneling event. 
These fluctuations can, of course also be induced by 
the tunneling current. However, in order to explain the 
observed width of 300 — 400 meV, single charges have to 
fluctuate within 2 nm of the cluster. We cannot exclude 



that enough defect states including chargeable cracks are 
close to each cluster. Thus, further studies are required 
in order to pinpoint the reason for the unexpectedly 
large FWHM. 



IV. CONCLUSION 

In conclusion, we have studied the electronic properties 
of MgO films grown on Mo(100) as well as of Co nano- 
islands on top of them. Using field emission resonances 
observed in scanning tunneling spectroscopy at room 
temperature, we found fluctuations of the work function 
within a closed MgO film up to A<£ = 0.6 eV. These fluc- 
tuations could be attributed to positively charged de- 
fect states (maybe F+ color centers) close to the MgO 
surface. On the Co islands, we observed Coulomb os- 
cillations at room temperature with energy distances in 
excellent agreement with orthodox theory. Analyzing 
the Coulomb gap we deduced a spatially varying back- 
ground charge Qo, which probably is also related to the 
charged color centers. The widths of the Coulomb peaks 
(300 — 400 meV) is considerably larger than the instru- 
mental resolution. We speculate that a strong coupling 
between the electrons in the Co islands and the phonons 
of the MgO induces this increased width. 
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